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We observe bright matter- wave solitons form during the collapse of 85 Rb condensates in a three- 
dimensional magnetic trap. The collapse is induced by using a Feshbach resonance to suddenly 
switch the atomic interactions from repulsive to attractive. Remnant condensates containing several 
times the critical number of atoms for the onset of instability are observed to survive the collapse. 
Under these conditions a highly robust configuration of solitons forms such that each soliton satisfies 
the condition for stability and neighboring solitons exhibit repulsive interactions. 
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Bose-Einstein condensates (BECs) with tunable inter- 
atomic interactions have been the subject of intense the- 
oretical and experimental interest in recent years. Many 
aspects of these systems are now well understood. How- 
ever, several questions remain pertaining to the response 
of the BEC to a sudden change in the nature of the inter- 
actions from repulsive to attractive. Resolving such ques- 
tions is critical if a complete understanding of these basic 
atomic systems is to be developed and if related, more 
complex phenomena, such as the formation of molecular 
BECs and superfluidity in degenerate Fermi systems, are 
to be thoroughly explained. 

The macroscopic wave function of an atomic BEC 
obeys a nonlinear Schrodinger equation, known as the 
Gross-Pitaevskii (GP) equation. The nonlinearity results 
from the interatomic interactions whose nature and mag- 
nitude are described by the s-wave scattering length (a). 
Homogeneous condensates in free space are inherently 
unstable when the interactions are attractive (a < 0). 
However, in a harmonic trapping potential, the zero- 
point kinetic energy can stabilize the condensate pro- 
vided the number of atoms in the BEC (No) is less than 
a critical value (N critical)- 



Nq < N cr itical — k 



O-ho 
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where is the mean harmonic oscillator length which 
characterizes the kinetic energy in the trap, and k is a 
dimensionless constant known as the stability coefficient. 
The exact value of k depends on the ratio of the trap 
frequencies and for our trap we have measured it to 
be k = 0.46(6) @. 

In the vicinity of a Feshbach resonance, a depends 
sensitively on the magnitude of an externally applied 
magnetic field 0, 0, Ej|, allowing the magnitude and 
sign of the atom-atom interactions to be tuned to any 



value. We have previously exploited such a resonance 
in 85 Rb to investigate the collapse of the condensate fol- 
lowing a sudden change in the nature of the interactions 
from repulsive to attractive [6j. The resulting dynam- 
ics of the imploding "Bosenova" were both fascinating 
and unexpected. Most notable was our observation of 
an anisotropic burst of atoms that explodes from the 
condensate during the early stages of collapse, leaving 
behind a highly excited remnant condensate that sur- 
vives for many seconds. Curiously, we frequently ob- 
served that the number of atoms in the remnant conden- 
sate, N remnant, was significantly greater than N 'critical- 
It wasn't understood why this remnant BEC didn't un- 
dergo further collapse until the number of atoms remain- 
ing was less than the value of N critical predicted by the 
GP equation. It was particularly puzzling that the sta- 
bility condition in Eq.^ accurately predicted the onset 
of the collapse, but did not seem consistent with number 
of atoms in the remnant. In the work presented here, 
we explain this result by showing that the remnant is 
composed of multiple solitons. 

Solitons are localized waves that propagate over long 
distances without change in shape or attenuation. The 
existence of solitonic solutions is a common feature of 
nonlinear wave equations. Solitons are therefore perva- 
sive in nature, appearing in many diverse physical sys- 
tems. BEC systems described by the nonlinear GP equa- 
tion can support both dark (local minima in the con- 
densate wave function) and bright (local maxima) soli- 
tons depending on whether the interactions are repul- 
sive (a > 0) or attractive (a < 0), respectively. Dark 
solitons have been observed in repulsive 87 Rb conden- 
sates 0, H| , while bright solitons have been created us- 
ing one-dimensional (ID) 7 Li condensates with attractive 
interactions 0, E3- Bright matter wave solitons form 
when the attractive atomic interactions exactly balance 



the wave packet dispersion. Gap bright solitons have 
also been observed using 87 Rb condensates confined in 
an optical lattice 0. In this letter we show that dur- 
ing the collapse of the condensate described in ref. 0, 
a robust configuration of multiple solitons forms in our 
three-dimensional (3-D) trap. Comparison with simu- 
lations shows that, remarkably, the solitons are created 
with set relative phases such that they repel each other 
and oscillate in the trapping potential for many seconds 
without degradation. Since the number of atoms in each 
soliton never exceeds N cr i t i ca i and the repulsive solitons 
never overlap, the condition for condensate stability is 
never violated. 

The experimental apparatus and the magnetic field 
ramps for probing the condensate collapse have been 
described in detail elsewhere 00]. Stable 85 Rb con- 
densates were produced within the positive scattering 
length region of a broad Feshbach resonance by radio- 
frequency evaporation in a magnetic trap. The frequen- 
cies of the cylindrically symmetric cigar-shaped trap were 
^radial = 17.3 Hz and Vaxiai = 6.8 Hz. Condensates 
of up to 15000 atoms were formed with a shot-to-shot 
standard deviation of less than 600 atoms. The non- 
condensate thermal fraction was less than 500 atoms. 
After the evaporation at 162.2 Gauss, the magnetic field 
was adiabatically increased to 165.5 G where a = +9ao 
(ao = 0.0529 nm is the Bohr radius). The collapse of 
the condensate was then induced by quickly switching 
the scattering length through zero to a variable negative 
value a co ii apse by increasing the magnetic field in 0.1ms. 
After a variable hold time t evo i ve at the negative scat- 
tering length, one of two imaging procedures was used 
to interrogate the remnant condensate. In the first, the 
scattering length was increased to a large positive value 
in order to use the large repulsive atomic interaction to 
expand the condensate prior to the standard rapid trap 
turnoff followed by absorption imaging sequence ;6] . This 
procedure ensured that the size of the condensate was 
well above the resolution limit of the imaging system 
(7/Ltm full- width at half- maximum (FWHM)), allowing 
the number of atoms to be measured reliably. In the sec- 
ond, the magnetic trap was turned off immediately and 
the condensate imaged with minimum delay (~ 2 ms), ef- 
fectively probing the size and shape of the remnant con- 
densate without any expansion. 

The number of condensate atoms surviving the col- 
lapse was investigated as a function of both a co ii a pse and 
the initial population of the condensate. The hold time 
at acoiiapse was sufficiently long (t evo i ve = 50 ms) to en- 
sure that all collapse-related atom loss had ceased |6|. 
The results are plotted in Fig.^ as the ratio of the rem- 
nant number to the critical number (calculated accord- 
ing to Eq.nfor each scattering length). Depending upon 
the conditions of the collapse, N remnant varied from be- 
low the critical number to much greater than the critical 
number (for a co u apse = — 100 ap, the remnant conden- 
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FIG. 1: Dependence of the number of condensate atoms sur- 
viving the collapse on the magnitude of the negative scatter- 
ing length at the collapse field (acoiiapse measured in units of 
the Bohr radius, ao = 0.0529 nm). The data are expressed as 
the ratio of the number of remnant condensate atoms to the 
critical number calculated according to Eq.Qfor each scatter- 
ing length. Results are shown for three different initial num- 
bers of condensate atoms; No = 15000 (•), 10000 (O) and 
5000 (a). The dashed line shows the critical number. The 
error bars represent the statistical spread in the data only. 



sate contained approximately 10 x the critical number 
of atoms). Interestingly, the fraction of atoms surviv- 
ing the collapse depends only on the scattering length 
and is independent of the initial number of condensate 
atoms (in the range we investigated) . This fraction varies 
smoothly from ~ 60% for a co llapse — ~5ao to ~ 30% for 
o-coiiapse — — 50ao. The lifetime of the remnant conden- 
sate at a co iiapse was several seconds and comparable to 
the lifetime of a condensate for small positive scattering 
lengths, even for N remnant > N criUca i. 

In a separate set of experiments, the evolution of the 
size and shape of the unexpanded remnant condensate 
was investigated. From the images, the axial and radial 
column density distributions were calculated and fitted 
with Gaussian profiles to obtain both the axial and ra- 
dial FWHM. The result of this analysis is shown for the 
(weaker) axial direction in Fig. [21(a). Any variation in the 
width of the remnant condensate in the radial direction 
was below the resolution limit of our imaging system. 
The remnant condensate is highly excited, with the vari- 
ation in the axial width being well fitted by sinusoidal 
oscillations at the two lowest collective mode frequencies 
(v\ = 2 X v axia i and v 2 = 2 x v radia i). The measured 
frequencies were v\ — 13.67(1) Hz and v 2 = 33.44(3) Hz. 
The amplitude of the v 2 oscillation was 29(7)% as large 
as that of the v\ oscillation and contained a significant 
contribution from the radial focus of the atom burst as- 
sociated with the collapse 0. There was no discernible 
damping during our 3 second observation time. 

Close inspection of the images, however, revealed that 
the observed variation in the width is not due to the ex- 
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FIG. 2: (color online) Observation of solitons oscillating in the 
magnetic trap following the collapse at a scattering length of 
o-collapse = —11-4 do of condensates initially containing ap- 
proximately 8000 atoms, (a) The evolution of the axial (hor- 
izontal) FWHM of the remnant condensate obtained from a 
single Gaussian fit to the images. Above the resolution limit 
of the imaging system, the remnant condensate is observed 
to separate into two solitons as shown in the images taken 
at (b) 210 ms, (c) 1140 ms and (d) 3110 ms. Each image is 
77 x 129 fim. The error bars represent the statistical spread 
in the data only. 



citation of collective modes. The remnant condensate is 
observed to separate into two (or more) distinct clouds 
which we associate with solitons (Fig.|S|(b-d)). The ob- 
served oscillation in the axial width of the remnant then 
results from the oscillation of two solitons along the ax- 
ial direction of the trap. As was observed for the soliton 
trains created in the 7 Li experiments [tj, neighboring soli- 
tons form with a relative phase, </>, that ensures that they 
interact repulsively (n/2 < <j) < 3tt/2) even though the 
atomic interactions are attractive [l^l H3|. In this case, 
the solitons do not pass through each other at the center 
of the trap, but rather rebound off one another. Conse- 
quently, as the solitons never fully overlap, the critical 
density for collapse is never reached and the solitons re- 
main stable. This conclusion is supported by the obser- 
vation that the soliton structures persist for over three 
seconds (Fig. 121(c)) with only a small degradation in the 
quality of the images due to the presence of more thermal 
atoms as the condensate melts. The solitons collide over 
40 times during this period. Several authors propose that 
bright matter wave solitons form due to a modulation in- 
stability of the condensate wave function 0,0,0. Our 
previous observations of local 'spikes' in the condensate 
density are consistent with this mechanism. 

The number of solitons that form depends both on the 
magnitude of a co ;/ apse and the initial number of conden- 
sate atoms. In Fig.|21we present images and cross sections 
showing two (Fig. [31(b)) and three (Fig.Q|(c)) solitons. 
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FIG. 3: (color online) Images and cross-sections of remnant 
condensates, (a) When the magnitude of a co ii a pse is suffi- 
ciently small a single remnant condensate containing less than 
the critical number is observed to survive the collapse. When 
the magnitude of a co iia P se is larger and/or larger initial con- 
densates are used, the remnant condensate is observed to split 
into a number of solitons determined by the conditions of the 
collapse (b-d). Each image is 77 x 129 /im. 



These images were acquired with t evo i ve = 200 ms, when 
the axial width was a maximum and the solitons were 
maximally separated. For comparison, Fig.El(a) shows 
a remnant condensate for which the number of remnant 
atoms was below the critical number. Under these con- 
ditions the remnant was never observed to separate into 
two or more solitons. We have observed up to six soli- 
tons present (Fig.|j3(d)). However, the quality of the 
images was significantly reduced when more than three 
solitons were present because the spacing between soli- 
tons became comparable to the imaging resolution limit. 
In Fig.EI(c) the three solitons contain approximately the 
same number of atoms. This was not always the case, 
however, and we often observed an unequal and (usually) 
symmetric population distribution across the solitons. 

In Fig.|H we examine in more detail the conditions 
under which multiple solitons form. For a fixed initial 
number of atoms in the condensate (No » 8000), the 
number of solitons created was investigated as a func- 
tion of a co ii apse . The images were again acquired with 
tevoive = 200 ms. In Fig.|H(a) the modal number of soli- 
tons is plotted. The data show a gradual increase in the 
number of solitons as the magnitude of a co ii apse increases. 
The ratio of the number of remnant condensate atoms to 
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FIG. 4: (a) The modal number of solitons (□) observed after 
200 ms of evolution at the collapse field as a function of the 
magnitude of the negative scattering (a co ;i apse ). The ratio 
of the number of remnant condensate atoms to the critical 
number (•) is also plotted on the right axis (note the two y- 
axes have the same scale). The initial condensate contained 
approximately 8000 atoms, (b) The ratio of the number of 
remnant condensate atoms to the critical number is plotted 
against the mean number of observed solitons to show that 
each soliton contains less than the critical number of atoms 
defined by Eq.^and is therefore expected to be stable. The 
error bars represent the statistical spread in the data only. 



the critical number for the same Nq is also plotted in 
Fig.@|(a). This ratio exhibits the same general trend as 
a function of a co u a pse as the observed number of soli- 
tons. Most importantly, it never exceeds the number of 
solitons. From this observation we deduce that the num- 
ber of atoms in each soliton is always less than or equal 
to the critical number. Each soliton is therefore stable 
against collapse, explaining the unexpected stability of 
remnant condensates for which the number of atoms is 
greater than the critical number. This point is further 
illustrated in Fig.0|(b) where the ratio of the number of 
remnant condensate atoms to the critical number is plot- 
ted against the observed mean number of solitons. All 
the points lie on or below the critical condition defined 
by Eq.Q] (dashed line). 

In contrast to the previous observations of bright soli- 
tons in one-dimensional confining potentials |9l ll0l| . these 
experiments were performed in a 3D trap with an aspect 



ratio of v r adiai/v a xiai ~ 2.5. Under such conditions, 3D 
bright solitons are predicted to be stable only if the at- 
tractive atomic interaction is smaller than a critical value 
|l4|. consistent with our observations. Simulations of our 
experiment based on a numerical integration of the 3D 
GP equation show that multiple solitons are only stable 
provided they interact repulsively Attractive inter- 
actions between the solitons would lead to energetic exci- 
tations of the soliton wave function during each 'collision' 
(in this case the solitons pass through one another), and 
ultimately a loss of particles from the BEC. Moreover, 
when the relative phase between solitons is fixed at 7r in 
the simulation, the maximum number of solitons that is 
predicted to be stable at a particular scattering length is 
in good agreement with our observations shown in Fig.0J 
The sudden change in the nature of the atomic inter- 
actions in a Bose-Einstein condensate from repulsive to 
attractive leads to a violent collapse process, in which a 
significant proportion of the initial condensate is ejected 
in a highly energetic burst. Yet, surprisingly, out of this 
violent explosion emerges a robust configuration of soli- 
tons that preserves a large fraction of the initial con- 
densate population. It is particularly striking that this 
behavior holds over a large range of attractive interaction 
strengths. 
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